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SPIFFI observations of the starburst SMM J14011+0252:i 
Already old, fat, and rich by z = 2.565 
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ABSTRACT 

Using the SPectrometer for Infrared Faint Field Imaging (SPIFFI) on the 
ESO VLT, we have obtained J, H, and K band integral field spectroscopy of the 
z = 2.565 luminous submillimeter galaxy SMMJ14011+0252. A global spectrum 
reveals the brighter of this spatially resolved system's two components as an in- 
tense starburst that is remarkably old, massive, and metal-rich for the early epoch 
at which it is observed. We see a strong Balmer break implying a > 100 Myr 
timescale for continuous star formation, as well as nebular emission line ratios 
implying a supersolar oxygen abundance on large spatial scales. Overall, the sys- 
tem is rapidly converting a large baryonic mass into stars over the course of only 
a few hundred Myr. Our study thus adds new arguments to the growing evidence 
that submillimeter galaxies are more massive than Lyman break galaxies, and 
more numerous at high redshift than predicted by current semi-analytic models 
of galaxy evolution. 
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1. Introduction 

Deep imaging in the optical and near-IR has made it possible to constrain both the evo- 
lution of the cosmic star formation rate density (e.g., Madau et al. 1996) and its time integral, 
the growth of the cosmic stellar mass density (e.g., Dickinson et al. 2003). Short-wavelength 
studies give an incomplete picture of these trends, however, since an important population 
of high-redshift galaxies is too obscured to be easily detected. These submiUimeter galaxies 
(SMGs; see Blain et al. 2002 and references therein) have large IR luminosities that are pre- 
dominantly generated by star formation rather than accretion (Barger et al. 2001; Almaini 
et al. 2003); as a result, they contribute substantially to cosmic star formation (Barger et 
al. 2000). Moreover, as the strikingly different appearances of the Hubble Deep Field at 
0.83 fim (Williams et al. 1996) and 850 fira (Hughes et al. 1998) exemplify, SMGs are rarer 
and forming stars much more intensely than typical optically selected systems. 

To understand the evolutionary state of SMGs, we have conducted a detailed rest-frame 
optical case study using SPIFFI, the SPectrometer for Infrared Faint Field Imaging (Tecza et 
al. 2000; Eisenhauer et al. 2000, 2003). Our target is SMM J14011+0252 (hereafter J14011), 
a. z = 2.565 galaxy lying behind the z = 0.25 cluster A1835 and the second SMG to have an 
optical counterpart identification (Barger et al. 1999) validated by CO interferometry (Frayer 
et al. 1999). Recent maps of J14011's molecular gas have led to divergent conclusions about 
its intrinsic size: Ivison et al. (2001) suggest that a large background source undergoes 
a factor A4 ~ 2.5 magnification by A1835 as a whole, while Downes & Solomon (2003) 
argue that a fortuitously located cluster member further boosts magnification of a small 
background source to ~ 25. We defer a full discussion of lensing to a future paper; here, 
we focus primarily on the lensing-independent conclusions of high age and metallicity that 
can be drawn from JMOll's global spectrum alone. Throughout the paper we assume a fiat 
0,A — 0.7 cosmology with Hq — 70kms~^ Mpc~^. 

2. Observations and data reduction 

In three SPIFFI "Guest Instrument" runs between February and April 2003, we observed 
J14011 in the J, H, and K bands on Kueyen (UT2 of the ESO VLT). Weather conditions 
were extremely stable within and between nights, and the near-IR seeing was better than 
0.6". SPIFFI obtains simultaneous spectra for all pixels in a 32 x 32 array; our data were 
obtained with the 0.25" pixel scale, and with spectral resolving power R = A/AA ~ 1500 
at 1.3 //m, ~ 2000 at 1.6 //m, and ~ 2400 at 2.2 //m. Individual exposure times in the three 
bands were 10, 5, and 10 minutes, respectively; we obtained an "off" frame at a sky position 
(offset by 20" — 25") for each "on" frame, and jittered the telescope by two-pixel offsets 
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between the "on" frames themselves. The total on-source integration times were 60 minutes 
in J, 95 minutes in H, and 340 minutes in K. 

We reduced the data within IRAF using modified versions of long-sht spectroscopy 
tools. After subtracting dark frames and fiat- fielding all data (object and sky frames), we 
masked out bad pixels and interpolated new values using all available three-dimensional 
information. Object and sky frames were then wavelength calibrated using arc lamp spectra 
and/or night sky emission lines. The "off" sky frames were then subtracted, and a linear fit 
to the background was subtracted across each spectral and spatial row of the data. At this 
point the data cubes for each on-off pair were reconstructed, using only integer pixel shifts 
and assuming that each spatial row of the field covered exactly 32 pixels; although this was 
not strictly correct, it only slightly reduced the spatial resolution and avoided loss of signal- 
to-noise due to resampling. A final background subtraction exploited the two-dimensional 
field of view: objects in the frame were masked out, and the mean of the remaining pixels 
in a given wavelength plane was subtracted from that plane. Using telescope offsets, we 
spatially aligned and then combined the data cubes, clipping any deviant pixels in each 
cell. Telluric correction was applied to the final merged cubes. Fluxes were calibrated by 
integrating SPIFFI observations of the UKIRT faint standard FS135 only over the Mauna 
Kea Observatory J, H, and K filter bandpasses; the reference photometry of FS135 itself 
has a scatter of less than 1% (Hawarden et al. 2001). 



3. Results 

Figure 1 shows J14011 in continuum-subtracted Ha and line-free K-band continuum. 
Within the bright eastern "Jl" component, the SPIFFI data reveal substantial spatial varia- 
tion in Ha equivalent width: nebular emission is strongest away from the central continuum 
peak. After using this peak to align all three data cubes, we extracted a spectrum of Jl 
from the 2" x 1.5" region shown in Figure 1. This box encloses about two thirds of the total 
continuum associated with Jl in the K-hand cube, and virtually all the emission we detect in 
J and H. Figure 2 shows the result and reveals that J14011 has a continuum break between 
the observed J and H bands. The lack of a similar break in blank-sky spectra extracted 
from elsewhere in the SPIFFI field proves that this cannot be an artifact of poor sky sub- 
traction; the J — H color of FS135, meanwhile, is accurate to 0.3% (Hawarden et al. 2001). 
The feature's most likely origin is the Balmer break of a > 100 Myr old stellar population 
at 2; ~ 2.565. Quantitatively, the strength of the break can be reproduced by a population 
synthesis model constructed with the STARS package (Sternberg 1998) for a 200 Myr old 
episode of continuous star formation with a 1 — 100 Mq Salpeter initial mass function (IMF), 
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solar metallicity, and an extinction Ay — 0.7. A crucial implication of this result is that 
the near-IR continuum emission must be dominated by the high-redshift background source, 
rather than by a foreground cluster galaxy. 

We have also used the spectrum in Figure 2 to measure fluxes and centroids for eleven 
nebular emission lines; taken together, they imply Jl has a systemic redshift zuii — 2.5652 ± 
0.0006, in perfect agreement with the zco derived from millimeter interferometry (Downes 
& Solomon 2003). H« emission from the fainter western "J2" component is blueshifted by 
Az = 0.0021 ± 0.0003 with respect to Jl, consistent with the offset measured by Ivison et 
al. (2000) from rest-UV spectra. The Jl line fluxes allow us to measure the i?23 estimator 
for its oxygen abundance (Pagel et al. 1979). After correcting for reddening with a Galactic 
extinction curve (Howarth 1983) and the color excess E{B - y) = 0.1811^;}^ implied by the 
Balmer decrement, and fixing 74959 = -^5007/3, we derive log(i?23) = 0.32 ± 0.16. Translation 
of this index to an oxygen abundance depends on ionization parameter, as quantified by 
log(032) = -0.24 ± 0.22. Using the fits of Kobulnicky, Kennicutt, & Pizagno (1999) to the 
calibrations of McGaugh (1991), we infer that 12 + log(0/H) = 8.96lo;?o- Relative to the 
measurement of Allende Prieto, Lambert, & Asplund (2001), this represents an abundance 
that is supersolar by 0.271q dex. A high [N II]/Hq; ratio means J14011 cannot fall on the 
lower-metallicity branch of the [0/H] vs. R23 relation. 



4. The high mass of J14011 

Although we reserve a full lensing model of J14011 for a future paper, our SPIFFI data 
already allow us to evaluate the small source/high Ai and large source/low Ai scenarios. As 
noted above, the strength of the apparent Balmer break in Jl leaves little room for continuum 
emission from an "extra" lensing galaxy in A1835. The newly confirmed velocity difference 
between Jl and J2, on top of Lya and GUI] 1909 A equivalent width differences in their 
pubhshed spectra (Barger et al. 1999; Ivison et al. 2000), indicates that these components are 
not multiple images of the same background source, contrary to the suggestion of Downes & 
Solomon (2003). Finally, SPIFFI spatially resolves Jl in an east-west direction; the observed 
0.95" FWHM of Ha emission exceeds the 0.65" east-west PSF dimension. This would not 
occur if Jl were formed from the merging of a triple image set. 

These direct arguments for a low lensing magnification are nicely complemented by the 
wealth of self-consistent evidence that J14011 as a whole has a large mass and age. Our 
estimate of Jl's high oxygen abundance is independent of Ai, provided there is minimal 
differential lensing from 3726 to 5007 A rest wavelength. For a closed-box model with solar 
yield, this metallicity corresponds to a gas/baryonic mass fraction of exp {—Z/Zq) ~ 0.16. 
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Prayer et al. (1999) measure a 8.2 x 10^° Al"^ K kms"^ pc^ CO(3-2) line luminosity for 
J14011. For gas that is dense and warm enough to produce detectable CO(7-6) emission, 
Downes & Solomon (2003) suggest that molecular gas mass and CO luminosity are related 
by Mgas/-^^co(3-2) ~ 0.8 M0 (Kkms~^ pc^)^-*^. Adopting --^ 5 as a fiducial value for the 
cluster-only lensing case, we infer a molecular gas mass Mgas ~ 1.3 X 10^° Mq, and therefore 
a total baryonic mass Mbary — 7.9 x 10^° Mq and stellar mass Mgtars = A/bary — Mg^s — 
6.6 x 10^° Mq. This last value (and with it, our assumption of moderate A4) is empirically 
supported by the mass-metallicity relation, as recently rederived by Tremonti et al. (2004) 
from a sample of ~ 51, 000 galaxies at 2; ~ 0.1. Their favored translation of R23 to metaUicity 
gives 12 + log(0/H) ~ 9.1 for Jl, at which the median stellar mass is ~ 6.2 x lO^'^ Mq. 

For a simple closed-box model, assumption of continuous star formation and a return 
fraction R — 0.2 uniquely sets the length of time for which a galaxy has been forming 
stars: isF = {exp (Z/Zq) - 1} (1 - R)~^ M^^ SFR~\ For J14011, ^sf will be independent of 
lensing because Mgas and SFR depend identically on Ai. We can estimate the star formation 
rate in J14011 from its ^850 = 12.3 A^-^mJy (Small et al. 2002); for the fiducial SED of 
Blain ct al. (1999) with Td = 40K, /3 = 1.5, and a = 2.2, Lir ~ 1.9 x 10^^ (^ggg/nijy) Lq = 
2.3x 10^'^ Al^^ Lq. Adjusting the conversion factor of Kennicutt (1998) to match a 1 — 100 Mq 
Salpeter IMF, we infer SFR ~ 1920 Mq yr~^ and tsF ^ 220 Myr. This age is consistent 
with the strength of the Balmer break seen in Figure 2, suggesting that the closed-box 
assumption made above is not unreasonable. J14011's powerful starburst thus appears to be 
converting most of its baryons into stars on a fairly rapid (few hundred Myr) timescale. 



5. The high masses of SMGs 

Its star formation history and metaUicity allow J14011 to join SMM J02399-0136- whose 
Mbary Gcuzel ct al. (2003) infer from its spatially resolved gas kinematics- as an SMG whose 
large baryonic mass is secure. Thanks to CO confirmation of Keck/LRIS redshifts obtained 
by Chapman et al. (2003b), this category promises to expand rapidly. Neri et al. (2003) have 
pubhshed CO data confirming the Keck redshifts for three new SMGs; the total set of five has 
{z) = 2.72 and a median CO linewidth of 420kms~^ (FWHM). A fruitful comparison can be 
made between this SMG sample and the well-studied Lyman break galaxy (LBG) population 
at z ~ 3. In the list of sixteen LBGs for which Pettini et al. (2001) report ionized gas 
kinematics, the median velocity width is 170 km s"-*^ (FWHM) at a slightly larger (z) = 3.12. 
Assuming that the (few kpc) size scales of the starbursts are comparable in LBGs and SMGs, 
the higher velocity widths imply SMGs' dynamical masses should also be larger by a typical 
factor of 6. If the two populations' dark matter fractions are comparable, we would expect a 
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median SMG baryonic mass exceeding the median 1.1 x 10^° M© stellar mass estimated for 
74 LBGs (Shapley et al. 2001)^ by the same factor, i.e., 6.6 x 10^° Although the strong 
clustering of LBGs (e.g., Giavalisco & Dickinson 2001) suggests they share their host halos 
with much more baryonic mass than is seen in their inner regions, SMGs appear already 
to have assembled measurably larger masses of baryons, or to have been more successful in 
retaining those baryons in the face of energetic winds. Complementary constraints on the 
star formation histories and metaUicities of more SMGs, like those obtained here for J14011, 
would further strengthen this case. 

In the hierarchical merging scenario for structiuc formation, massive galaxies are as- 
sembled in greater numbers with the passage of time. As several authors have pointed out, 
the observed number densities of high-redshift galaxies with large stellar masses exceed the 
predictions of semi-analytic models of galaxy evolution within this framework (Cimatti et 
al. 2002; Daddi et al. 2004; Saracco et al. 2004). As in Genzel et al. (2003), we can also 
use the number densities of SMGs to apply the same "baryonic mass assembly test" to the 
models. Small et al. (2002) report the surface density of SMGs with 8350 > 5 mJy (the 
threshold above which SMG optical redshifts 1.8 < z < 3.5 have been confirmed with CO 
interferometry) to be 88OI33Q deg^^. The Chapman et al. (2003b) sample includes two (of 
nine) sources above this threshold lying at z < 1.8, and excludes the ~ 35 ± 5% of bright 
SMGs that lack radio counterparts (Ivison et al. 2002; Chapman et al. 2003a). Two of the 
five SMGs with CO detections have > 510kms~^ FWHM line widths that would correspond 
to Mbaxy > 10^^ Mq by the above rescahng of LEG stellar masses. This leaves an estimated 
surface density of lSO't\l^ deg~^ for SMGs with baryonic masses > 10^^ Mq in the range 
1.8 < 2; < 3.5, and a comoving number density of S.QlgJ x 10~^Mpc~^. 

We must still correct this density upward to account for sources that are comparably 
massive but no longer luminous enough to have been detected at 850 /im. Adding ^sf and 
the gas exhaustion timescale Mgas SFR~^ for J14011 implies that the total duration of its 
current burst will be of order 250 Myr. Relative to the 1.8 Gyr elapsed over the range 
1.8 < z < 3.5, this would imply a correction factor of 7 for the population as a whole. 
Figure 3 plots this prediction, together with the comoving number densities of galaxies with 
stellar masses above the 10^^ Mq threshold at 2; ~ (Cole et al. 2001) and z 1 (Drory et al. 
2004). Also plotted are the predictions of semi-analytic models, one of which is taken from 
Kauffmann et al. (1999) and the other an update of Baugh et al. (2003). Since all theoretical 
and observational values assume the same cosmology and a consistent (either Miller-Scalo 
or 1 — 100 Mq Salpeter) IMF, we are secure in concluding that the models underpredict 



^This median LBG mass adjusts the Shapley et al. (2001) population synthesis results for a 1 — 125 
Salpeter IMF and Hq = 70. 
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the number densities of massive galaxies at 2; ~ 2.7. Couching Figure 3 in terms of mass 
emphasizes that the observed surface densities of SMGs cannot be explained using a very 
fiat IMF alone. We also note that the models contain sufficiently many dark halos of total 
mass > 10-*^^ (Qi/i/^h) Mq to account for the observed number densities of massive SMGs, 
provided their baryons can be rapidly and efficiently assembled into galaxies. 

Correcting the observed number density of SMGs to account for their less dust-luminous 
descendants raises the question of what types of galaxies SMGs can plausibly evolve into. 
Although LBGs are clearly excluded due to their smaller masses, an intriguing alternative is 
the newly identified population whose red J — K colors can stem from a strong Balmcr break 
at 2; ~ 2.5 (Franx et al. 2003). These galaxies appear to be strongly clustered (Daddi et al. 
2003) and- assuming a uniform distribution in the redshift range 2 < z < 3.5 (van Dokkum 
et al. 2003)- have a comoving number density ~ 1.8 x 10~^Mpc~^. Relatively few sources 
with J — X > 2.3 also have 1.2 mm flux densities > 3mJy (Dannerbauer et al. 2004), 
indicating little direct overlap with the bright end of the SMG population (see also Frayer et 
al. 2004), although fainter (sub) millimeter counterparts may still be present (e.g., Wehner, 
Barger, & Kneib 2002). With three examples in the HDF-S having (for a Miller-Scalo IMF) 
stellar masses (0.6 — 1.4) x 10^^ Af© (Saracco et al. 2004), it would seem plausible that the 
more luminous objects with red J — K colors at this epoch could have passed through an 
SMG phase at higher redshift. 

We thank M. Horrobin, C. Iscrlohc, A. Schegerer, and J. Schreiber for their hard work 
on SPIFFI; A. Sternberg for assistance with STARS; C. Baugh, N. Drory, G. Kauffmann, 
and C. Tremonti for useful discussions; and an anonymous referee for helpful comments. 
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Fig. 1. — Left: continuum-subtracted Ha image of J14011. Right: line-free X-band contin- 
uum image of J 140 11. The field of view is 8" x 8"; the aperture used to extract the composite 
spectrum for Jl in Figure 2 is overlaid. 
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Fig. 2. — J, /J, and K spectra of J14011 extracted from a 2" x 1.5" aperture. Superposed is 
the model STARS spectrum for a 200 Myr old continuous star formation episode, assuming 
solar metaUicity, a 1 — 100 Mq Salpeter IMF, and extinction Ay — 0.7. 
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Fig. 3. — Comoving number densities of galaxies with baryonic masses > 10^^ Mq as a 
function of redshift. Triangle and open squares show densities of massive stellar systems at 
2 = (Cole et al. 2001) and 2; ~ 1 (Drory et al. 2004); circle shows density for massive 
SMGs at 2; ~ 2.7, with a factor of 7 correction for burst lifetime (see §5). Solid curves show 
the predictions of semi-analytic modelling by the "Munich" group (Kauffmann et al. 1999) 
and the "Durham" group (an updated version of Baugh et al. 2003); dashed curves show 
the corresponding number densities of halos with available baryonic masses > 10^^ Mq. The 
two models use the same halo simulations but assume different Qb- 



